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Abstract 

Discrete nonlinear Schrodinger equation (DNLS) describes a chain of oscillators 
with nearest neighbor interactions and a specific nonlinear term. We consider its 
modification with long-range interaction through a potential proportional to l/l l+a 
with fractional a < 2 and I as a distance between oscillators. This model is called 
aDNLS. It exhibits competition between the nonlinearity and a level of correlation 
between interacting far-distanced oscillators, that is defined by the value of a. We 
consider transition to chaos in this system as a function of a and nonlinearity. It 
is shown that decreasing of a with respect to nonlinearity stabilize the system. 
Connection of the model to the fractional genezalization of the NLS (called FNLS) 
in the long-wave approximation is also discussed and some of the results obtained 
for aDNLS can be correspondingly extended to the FNLS. 
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1 Introduction 

Nonlinear Schrodinger equation (NLS) is a paradigmatic equation that describes a slowly varying 
enveloping process in the nonlinear dispersive media. Applications of the NLS has been found 
in almost all important areas of physics: nonlinear optics, plasma physics, hydrodynamics, 
condenced matter physics, biology and others. The literature on NLS is extensive and the 
reviews ^ |2J El HI 00 E] can provide a strong impression on the importance of the subject. The 
NLS per se is an integrable system while its different types of perturbations, more related to the 
practical needs, include a broad spectrum of solutions from solitons and breathers to chaos and 
spatio-temporal turbulence. 
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Between different types of perturbations one can single out two the most interesting classes: 
external time-space dependent perturbation, and discretization of the NLS presented by a kind 
of space-difference equation, called discrete NLS (DNLS). It was shown in a set of publica- 
tions (see (7J |HJ and references therein) that specific hyperbolic structure of the NLS phase 
space leads to the transition from a soliton type dynamics to different kind of chaos, including 
spatio-temporal turbulence, under even small perturbation. The fact that discretization induces 
chaos is fairly well known and, particularly, for the DNLS it was studied in |10| 1121 113| I14j 
and |151 I16j . Destruction of solitons, breathers, and wave trains is similar to what occurs from 
external perturbations. A general physical mechanism of the onset of chaos induced by a dis- 
cretization is known: transition from a differential equation to the difference one is equivalent 
to the appearance of a high frequency (in time or in space) periodic perturbation |17j . 

One can also consider DNLS as a separate problem that describes a chain of coupled many- 
particles (oscillators) with local or nonlocal interaction. Such a system presents a specific interest 
for studying transition to chaos, turbulence, and statistical equilibrium in many-body problem. 
A new attraction at that point is the long-range interaction (LRI) that was introduced in an 
exponential Kac-Baker form in |18| I19j and later in a power Lennard- Jones form in [201 • ^ n 
another version of the latter case the interaction between oscillators located at the positions 
(n,m), n ^ m, is proportional to l/\n — m\ 1+a . The corresponding DNLS will be called aDNLS. 
Studying of such systems has multiple interest: transition to chaos and turbulence in the presence 
of LRI |21| I22j , sinchronization in systems with many particles [22] , controling chaos 24 , and 
different manipulations with physical objects in optics PJ EH] and condenced matter |26[ I27j . 
All these physical features, being described by aDNLS, are functions of a. 

For the number of particles N = 4 the DNLS equation can be solved exactly. For N > 4 
the DNLS is not integrable and chaotic solutions are possible. It was shown in |131 1151 1281 
I29 [ \'AU \ l3~T] that there are two main mechanisms responsible for chaotization of the DNLS. One 
of the mechanisms prevails for symmetric initial conditions and another for asymmetric ones. 
For symmetric initial conditions chaos was shown to emerge from the proximity to homoclinic 
orbits. For asymmetric initial conditions this mechanims do not play significant role, instead, 
perturbations induced by discretization and round-off errors cause random flipping of wave's 
direction of motion. 

The primary goal of this paper is to study the transition to chaos in aDNLS depending on 
a parameter < a < 2 that is responsible for the appearance of a power-like tail in solution, 
i.e. for the derealization of modes. Parameter a has a simple physical meaning: it describes a 
level of collective coupling of particles. Particularly, for a = oo we have only nearest neighbor 
interaction and for a = — 1 we have mean field type model. 

Recenly, different properties of the nonlocal DNLS were studied in [321 EHH EH OH! ■ ^ was 
found that some properties of solutions depend on the interaction exponent a in the LRI case or 
on the exponent (3 for the Kac-Baker interactions. Namely, for the power-law interaction it was 
shown that for a less than some critical value a cr , there is an interval of bistability where three 
posible stationary states exist at each value of some excitation number M that characterize 
a level of nonlinearity. The first type of solution is a continuum-like mode, the second is an 
intrinsically localized breathing state and the third one is stationary itermediate state. The 
long-distance behavior of the intrinsically localized states depends on a. For a > 2 their tails 
are exponential, while for < a < 2 the tails are algebraic. By changing an interaction constant 
J, a stable solution may become unstable. A small symmetric force applied to the system can 
also trigger transigions from one type of stable solutions to another [33]. The continuum limit 
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of this model was shown to be a nonlocal NLS equation [32]. Another version of the nonlocal 
(integro-differential) NLS equation was proposed in [35]. Unlike the usual NLS equation, this 
nonlocal NLS equation has stationary solutions only in a finite interval of excitation numbers 
[0,M max ]. 

In a similar way, as dynamics of a chain of coupled oscillators can be reduced to the wave 
equation in the long wave-length limit k — * 0, the chain of aDNLS can be reduced to the 
fractional generalization of NLS equation or Ginzburg-Landau equation (FGL) [HZ1 EE1 EES 1221 
EQ1ET]. It was shown in jHHIlSj that mapping of the aDNLS equation to fractional NLS (FNLS), 
or similarly to FGL, can be realized by some transform operator. In all fractional equations of 
this type second coordinate derivative is replaced by the fractional Riesz derivative 02] of order 
a. The corresponding comparison of solutions of discrete chain of oscillators with the sine- 
potential and fractional generalization of the sine-Gordon equation was perfomed in |27[ I43| . 
The latter results of [13] confirm a path to the dual features of the systems with LRI of order a 
and the systems described by a corresponding equation with fractional derivative I2.3| . 

In this paper we provide a detailed study of the transition to chaos in aDNLS depending 
on the parameter < a < 2. It is well known that fractional values of a appears in numerous 
complex systems such as spin-interacting systems [II], adatoms |26| . colloids |45j. chemical 
surfaces [Hj, quantum field theory |I3 HB1 HH] , etc. In correspondence to the results [39 ] 123 ] HT)] . 
the obtained properties of the aDNLS can be extended to the FNLS. It will be shown that 
onset of chaos follows as a result of competition between the nonlinearity level and the level of 
coherency of the chain of oscillators defined by the value of a. 



2 Basic equations 

The continuous Nonlinear Schrodinger (NLS) equation defined on the finite interval [-L/2, L/2] 
with the periodic boundary conditions ip(x + L,t) = tp(x,t) can be written as 

*- +7 |^+_ = , (1) 

where 7 is a constant (7 = 1 corresponds to the focusing nonlinearity). In turn, the discrete 
nonlinear Schrodinger equation (DNLS) which describes a lattice of N anharmonic oscillators 
with nearest-neighbors interaction is defined as 

i^ + 7l^n|Vn + e(V'n+i - 2^„ + Vn-i) = 0, (n = l,...,N), (2) 

where ip n +N = ipn, Vn is the periodic boundary condition. The quantity ip n = ip n (t) is the 
complex amplitude of the oscillator at site n. With e = l/(Ax) 2 , Eq. (J2J) is seen as a standard 
finite difference approximation to Eq. Here Ax = L/N is a distance between oscillators. 
The Hamiltonian, H, and the excitation number (or norm), M, 

N N 

H = J2(e\i>n+l-Tpn\ 2 - \\^nf), M = £ |^| 2 , (3) 
n=l n=l 

are the conserved quantities. 
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The model which we study in the following is described by the Hamiltonian 

(4) 



1 N ( N \ 

n=l \ m=l I 



with the periodic boundary conditions: ip n +N = ipn- The coupling function is defined by 

Jn-m = i J Mir , (n + m), (5) 

\n — m\ L+a 

where J is a coupling constant and a is an exponent which, depending on the physical situation, 
can be integer or fractional. The standard nearest-neighbor DNLS equation is recovered in the 
limit a —* oo. From the Hamiltonian Eq. Q one obtain the equation of motion 

i^ = dH/thfc, (6) 
dip N 

i—^- +7l^n| 2 V'n+ ^2 J n-m(lpn ~ Ipm) = 0, (n = 1, N). (7) 
m=l 

Consider the infinite chain of equidistant oscillators (N — > oo in Eq. Q). The Fourier 
transform of ip n is given as 

oo 

4> n (k,t) = Vn(i) exp(-«fcc n ) = ^A{^n(t)}, (8) 

n=— oo 

where x n = nAx is a coordinate of the n-th oscillator, and Ax = 2ir/K is a distance between 
oscillators. Here we have treated A; as a continuous variable. The inverse transform is defined 
as 

1 f +K ' 2 . « 
V'n(i) = — I dk ip n (k,t) exp(ikx n ) = T A (9) 

-JC/2 



A" 



Transition to the continuous limit Ax — > (iT — > oo) can be obtained by transforming ip n (t) 
ip(nAx,t) = tp(x n ,t) — ► ip(x,t). Then, changing sums into integrals, Eqs. ©, © become 



$(M) = / ^x e~ ikx ip(x,t) = F{ip(x,t)}, (10) 

J — oo 
1 f + oo 

^(x,t) = — dk e ikx ip(k,t) = F~ x {i;{k,t)}, (11) 

where 

ip{k,t) = £ip(k,t), ip{x,t)=£ip n {t)=£ip(x n ,t), (12) 

and £ denotes the limit Ax — > 0. Operation T = J r ~ 1 £ can be called a transform operator 
(or transform map), since it performs a tranfrormation of a discrete model of coupled oscillators 
to the continuous media. For < a < 2 application of T leads to the fractional NLS. For more 
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details, definitions and applications of the transform operator to different systems see [2311401 IH?] . 
In a brief form this approximation is as follows |391 I23| . 

First we perform the Fourier tranformation of Eq. (J7J) for infinite chain of oscillators (N —* oo) 

d$(k,t) 

where 



+ 7 FAiWnfM + J (Ja(k) - J a (0)) ij(k, t) = 0, (13) 



I" 00 g—iknAx g-ifcnAx _|_ gjfcnAa 



Uk) = £ = £ ^ = ^ +a (e lkAx ) + L il+a (e~ lkA % (14) 

n=— oo ' ' n=l 

and Lzi +a (z) is a polylogarithm function JSUJ with series representation 

J a {k) = a a |Ax| Q |fc| Q + 2 V C(1 X" \T ^ (Ax) 2n (-k 2 ) n , \k\<l, a + 0,1,2,3... . (15) 

Here J Q (0) = 2£(1 + a), £ is the Riemann zeta-function and 

a a = 2 T(-a) cos (^) . (16) 

For a = 2, J a (fc) reduces to the Clausen function C^ik) [ST] . 
Substitution of Eqs. (|T5|) and (fl6|) into (|T3|) gives 



.dip{k,t) 



m + J ^k,t) (^a a \Ax\ a \k\ a + 2 £ C( " 1L. 2n ) (Ax) 2 "(-fc 2 )" - J Q (0)^ + (17) 

+7^A{|V>n|Vn} = 0. 

Note that J Q (0) exactly cancels the constant which is the first term of the sum in Eq. (|17|) . In 
the limit k -> Eq. dTTJ) yields 

+ j f ^ {k) ^ t) + ^ ^n^^ = 0) (18) 



where J = J\Ax\ min ^ and 



f fA A = JWi - |Ax| 2 - a C(a - 1)A: 2 , 0<a<2, (a ^ 1) 
"■ A ■' S \Ax\ a - 2 a a \k\ Q -((a-l)k 2 , 2<a<4, (a / 3). 



For < a < 2 the operator T aj A(k) is defined up to 0(k 2 ) and for 2 < a < 4 up to 0(|fc| a ). 
Eq. Q19|) has a crossover scale for 

fco = |a a /C(a-l)| 1/{2 ^ ) |Ax|- 1 . (20) 

From Eq. ^ it follows that T aA {k) ~ A; 2 for a > 2, A; <C ko and nontrivial expression 
^a,A(A:) ~ appears only for a < 2, /c <C /cq. This crossover was considered also in [3*91 125] . 
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Performing the transition to the limit k <C &o (or more precisely kAx <C koAx) and applying 
the inverse Fourier transform to (|18|) . we obtain 

d 

i—i(j(x,t) + JT a {x)i;(x,t)+ 7 \i;(x,t)\ 2 ij(x,t) = a + 0,1,2,..., (21) 



where 



i - \-a a -Jr^, 0<a<2, (a^l) , s 

r Q (x) = ^- 1 {r Q (fc)} = r/^; 2 n ' / \;/ (22) 



r a (fc) 



C(a-l)^, a>2, (a/ 2,3,4,...); 

'a a |/c| a , 0<a<2, (a^l) 

-((a-l)k 2 , a > 2, (a / 2,3,4,...). 



Here, we have used the connection between the Riesz fractional derivative and its Fourier trans- 
form 

— -%r * (23) 

Properties of the Riesz fractioanl derivative can be found in |42 | I52 | I53* | I54j . 

In the following section we will consider simulations of the finite chain of oscillators on a 
finite interval (— L/2,L/2) (n = 1, —,N, where N is even) described by the Hamiltonian 

ff = E JE h^zM .2^ (24) 

n=l \ 1=1 J 

and equation of motion 

N/2-1 

2 — +7|V>n| V>n + J 2^ =0 ' (« = l>->^)- I 25 ) 

1=1 

The form of the Hamiltonian and the equation of motion account for the periodic boundary 
conditions vp n +N = ipn- To avoid double counting of interactions we introduce a cutoff for the 
interaction range at N/2 — 1. 

All initial conditions and parameters are choosen to make it possible to use the transform 
operator T. This means that for not too large time the obtained results can be also applied to 
the FDNLS equation (J2U) 

d d a ib(x t) 

i-^(x,t)-G > %b(x J t)+ 1 \iP(x,t)\ 2 iP(x,t)=0, 0<a<2, a ^ 1, (26) 



G = 2 J \Ax\ a F(-a) cos(7ra/2). (27) 

It is also worthwile to mention that for a — > 2 we satisfy the conditions of the continuous limit 
approach since ko — > oo, that means the corresponding results should be close to the solution 
of NLS. Thus, considering a far from a = 2, one can compare the solutions of aDNLS with 
the solutions of DNLS with the nearest neighbor interaction. Such results provide a role of the 
long-range interaction for the chain dynamics. 



6 



3 Numerical results 

In this section numerical results obtained from solutions of the equations of motion Eq. (|25[) on 
the finite interval (— L/2,L/2) are summarized. Parameter J was normalized to J/ Jo, where 

N/2-1 

J o = E (28) 

n=l 1 1 

For all sets of parameters we integrate the equations of motion Eq. I|25jl up to the time T = 100. 
Numerical solutions were stored at each t„ = qT/Q, q = 0, .., Q — 1, where Q = 10 3 . The number 
of oscillators was N = 32 in all simulations. 

To visualise the numerical solution we plot the surface \ip(x, t)\ 2 and the "phase portrait" of 
the central oscillator (n = 0, x n = 0) formed by the variables: 

A(t) = |^(0, t) | 2 , A t = dA/dt. (29) 

Also we plot lm(ip(0,t)) vs. Re(ip(0,t)) and the phase difference of two nearby trajectories 
V>(0,t), ^(Ax,t) 

df = tan" 1 ( J Re(V'(0,t))//m(V'(0,t))) - tan' 1 (Rety (Ax, t))/Im(ip(Ax, t))), (30) 

where Ax = L/N. We also calculate the discrete Fourier transform of the sequence '0(0, tq) , 
which is defined 

1 

$n(wj) = qY1 ^(tq) exp(-« Wj t g ), (31) 

^ g=0 

Q-l 

i=o 

where the wavenumber is Wj = 2-Kj/Q, j = 0, ...,Q — 1. The corresponding power spectrum S 
of the sequence ip n {tq) (q = 0, Q — 1) is given by 

5,- = S^) = |^ n K-)| 2 - (33) 



In our simulation we use two types of initial conditions: symmetric and asymmetric types. 
Symmetric initial conditions are defined as 



V>(z,0) 



{2ir 
a + ocos I — x 



(34) 



with constants a = 0.5, b = 0.1 and L = 2^/2n. This form of symmetric initial conditions was 
also used in |15| l28j to study chaos in discretizations of standard NLS equation. For asymmetric 
initial conditions we use the following expression 



i>(x,0) 



l + b 



/2vr 



e"~cos I —x | + e tu- sin 



2vr 



(35) 



where a = 1, b = 0.2, c = 0.9, d = 2.03 and L = 2\f2n. The form of the initial condition is 
taken the same as in |13j . 
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Figure 1: Time evolution of the central oscillator. The values of parameters are a = 1.11, 
J/ Jo = 0.7 and M = 12.5. The initial condition is given by Eq. (34). 
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Figure 2: Time evolution of the central oscillator. The values of parameters are a = 1.11, 
J I Jq = 0.7 and M = 14.28. The initial condition is given by Eq. (34). 
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Figure 3: Time evolution of the central oscillator. The values of parameters are a = 1.11, 
J I Jq = 0.7 and M = 14.92. The initial condition is given by Eq. (34). 
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Figure 5: Surfaces \ip n (t)\ 2 for the same values of parameters as in Figs. 1-4: a = 1.11, 
J/J = 0.7 and M = 12.5 (top left), M = 14.28 (top right), M = 14.92 (bottom left) and 
M = 16.6 (bottom right). 
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Our main goal is to compare solutions of Eq. (|25|) for different values of a £ (1,2) and 
consider transition to chaotic dynamics of the chain of oscillators as a function of a. The larger 
is M, the stronger is nonlinearity. The larger is a, the smaller is LRI. The main differences in 
the physical properties of the oscillators dynamics are defined by a competition between a and 
M. In all simulations we fix J/Jq = 0.7 as the value close to the one considered in |32| . All 
plots will show the following properties of the oscillators: plane (dA/dt, A) shows projection of 
the trajectory of the central oscillator in phase space (see definition in (|29)0 : plane (Im z, Re z) 
shows projection of the complex amplitude z = ip(0,t) of the central oscillator as a function of 
time; phase difference with the adjacent oscillator to the central one (see Eq. (f30|) ): spectrum of 
time oscillations of ip(0,t) (see definition in Eq. (|3*3*|0 : and surfaces \ip n (t)\ 2 vs t and n. 

First four Figures 1-4 aim to show different regimes of the chain of oscillators behavior for 
a = 1.11. For small values of M solutions are quasi-periodic in time with only few modes in 
the Fourier spectrum. This behaviour for M = 12.5 is shown on Fig. 1. The plane (dA/dt, A), 
the plane (Im z, Re z) and the phase difference plane demonstrate quasi-periodic behavior. As 
we increase the value of M, the spectrum is broadening. Figure 2 demonstrates the behaviour 
of the system with a = 1.1, J/ Jo = 0.7 and M = 14.28. The quasi-periodic structure of the 
plane (dA/dt, A), the plane (Im z,Re z) and the phase difference plane get more complex. 
More and more Fourier modes appear in the spectrum. This is even more pronounced for 
the case M = 14.92 depicted in Fig. 3, which can be considered as a begining of chaos. The 
phase difference of two nearby oscillators shows two 'flips' to it and — ir which indicates phase 
decoherence and transition to chaos. In Fig. 4 for M = 16.6 the phase difference of two nearby 
oscillators has many 'flips' to ir and —it and the Fourier spectrum of ip(0, t) becomes broad what 
is typical for chaotic dynamics. Surfaces \ip(x,t)\ 2 for the cases of Figs. 1-4 are shown in Fig. 5. 
In addition to the chaotization of the solution which was described above, it is seen that with 
increasing values of M, surfaces become more localized around x = 0, oscillations in the wings 
of the solutions gradually decrease and the aplitude of the solution is increased. This can be 
explained by the growth of the nonlinear term with the growth of M. The role of a nonlinear 
coupling becomes more important than coherent connection of oscillators due to the LRI. 

Increasing of a leads to appearence of chaotic dynamics for smaller M, without significant 
changes of the diagrams and planes shown in Figs. 1-5. This is demonstrated for a = 1.91 in 
Fig. 6. For a = 1.51 chaos starts at approximately M = 14.28 and for a = 1.91 at M = 12.5. 
Dynamics for a < 1 is approximately similar to the dynamics for 1 < a < 2 and transition to 
chaos for a = 0.73 occurs at M ~ 18. The growth of a leads to the increasing of the energy of 
oscillations in the tails of the solutions. Transition to the strongly developed chaos is not too 
sharp in time. 

All these results could be compared to the standard DNLS equation with the nearest- 
neighbor interaction. Solution of this equation with the symmetric initial condition Eq. (|34[) 
is shown in Fig. 7 for M = 12.5. There are several main distinctions of this case from the 
transition to chaos when a < 2: (i) the symmetry of the solution breaks down for larger M and 
sharply in time; (ii) the spatial chaos occurs visually earlier than in the case of a < 2, and this 
indicates the role of LRI comparing to the standard case of the nearest-neighbor interaction; (iii) 
another important difference of the onset of chaos in aDNLS and in DNLS can be deduced from 
the (Im z, Re z) plots; in Fig. 7 (top right panel) for DNLS, trajectory fills space more-or-less 
uniformly what is typical for Hamiltonian chaos, while in Figs. 3, 4 for aDNLS trajectories 
looks like in the case of stochastic attractors what is more natural for a < 2 |56| . There exists 
an inner part of the diagram that is avoided by the trajectories, at least for the observed time. 
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Figure 7: Solution of the standard DNLS equation with M = 12.5 and J/ Jo = 0.7. The 
initial condition given by Eq. (34). 
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Figure 8: Time evolution of the system of coupled oscillators with LRI and asymmetric 
initial conditions. The values of parameters are a = 1.11, M = 22.22, J/ Jo = 1. The 
initial condition is given by Eq. (35). 
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The last considered case is the asymmetric initial condition given by Eq. (|35|) . Numerical 
solution of the equations of motion reveals a difference in this case compared to the symmetric 
initial condition case. For some values of the excitation number M the numerical solution starts 
to move in the left or right dirrection. This dirrection can also change randomly in time. Figure 
8 shows an example of such behavior for a = 1.11, J/ Jo = 1 and M = 22.22. Note, that the 
power spectrum in this case is also broad. 

4 Conclusion 

One of the main feature of the considered aDNLS model is implementation of a new, additional 
to the standard DNLS, parameter a that in the continuous limit implies the fractional dynamics 
described by the FNLS. From another point of view, a is responsible for stong correlations 
between distant oscillators, i.e. a long-range interaction is introduced through the parameter 
a. This feature of the aDNLS brings a new physics with a new control parameter. The role 
of the LRI was known before for collective phenomena in complex medium such as chemical or 
biological set of objects [S7j, phase transition in one-dimentional systems |58| . synchronization 
[2*H] . regularization in quantum field theory [IB]- Our detailed analysis helps to understand some 
specific properties of destabilization and onset of chaos in aDNLS with a < 2. Similar analysis 
can be performed for other models with LRI. An important part of our analysis is utilization 
of the possibility to transfrom the behavior of discrete chain of interaction objects into the 
continuous medium equation with the fractional derivatives. This formal procedure raise the 
question of the discrete-continuous equivalence up to a new level where the appearance of an 
additional parameter a increases the difficulty the answer the equation. 
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